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Many physically motivated extensions to general relativity (GR) predict sig- 
nificant deviations in the properties of spacetime surrounding massive neu- 
tron stars. We report the measurement of a 2.01+0.04 solar mass (Mq) pul- 
sar in a 2.46-hr orbit with a 0.172+0.003 M^ white dwarf. The high pulsar 
mass and the compact orbit make this system a sensitive laboratory of a pre- 
viously untested strong-field gravity regime. Thus far, the observed orbital 
decay agrees with GR, supporting its validity even for the extreme conditions 
present in the system. The resulting constraints on deviations support the use 
of GR-based templates for ground-based gravitational wave detectors. Addi- 
tionally, the system strengthens recent constraints on the properties of dense 
matter and provides insight to binary stellar astrophysics and pulsar recycling. 

Neutron stars (NSs) with masses above 1.8 Mq manifested as radio pulsars are valuable 
probes of fundamental physics in extreme conditions unique in the observable Universe and 
inaccessible to terrestrial experiments. Their high masses are directly linked to the equation- 
of-state (EOS) of matter at supra-nuclear densities (7, 2) and constrain the lower mass limit 
for production of astrophysical black holes (BHs). Furthermore, they possess extreme internal 
gravitational fields which result in gravitational binding energies substantially higher than those 
found in more common, 1.4Mo NSs. Modifications to GR, often motivated by the desire for 
a unified model of the four fundamental forces, can generally imprint measurable signatures in 
gravitational waves (GWs) radiated by systems containing such objects, even if deviations from 
GR vanish in the Solar System and in less massive NSs (3-5). 

However, the most massive NSs known today reside in long-period binaries or other systems 
unsuitable for GW radiation tests. Identifying a massive NS in a compact, relativistic binary 
is thus of key importance for understanding gravity-matter coupling under extreme conditions. 
Furthermore, the existence of a massive NS in a relativistic orbit can also be used to test current 
knowledge of close binary evolution. 

Results 

PSR J0348+0432 & optical observations of its companion PSR J0348+0432, a pulsar spin- 
ning at 39 ms in a 2.46-hr orbit with a low-mass companion, was detected by a recent sur- 
vey (6, 7) conducted with the Robert C. Byrd Green Bank Telescope (GBT). Initial timing ob- 
servations of the binary yielded an accurate astrometric position, which allowed us to identify 
its optical counterpart in the Sloan Digital Sky Survey (SDSS) archive (8). The colors and flux 
of the counterpart are consistent with a low-mass white dwarf (WD) with a helium core at a dis- 
tance of J ~ 2.1 kpc. Its relatively high apparent brightness (g' = 20.71 + 0.03 mag) allowed us 
to resolve its spectrum using the Apache Point Optical Telescope. These observations revealed 
deep Hydrogen lines, typical of low-mass WDs, confirming our preliminary identification. The 
radial velocities of the WD mirrored that of PSR J0348+0432, also verifying that the two stars 
are gravitationally bound. 



In December 201 1 we obtained phase-resolved spectra of the optical counterpart using the 
F0RS2 spectrograph of the Very Large Telescope (VLT). For each spectrum, we measured the 
radial velocity which we then folded modulo the system's orbital period. Our orbital fit to the 
velocities constrains the semi-amplitude of their modulation to be ^wd = 351 ±4 km s"^ (Fig. 1; 
see also Materials & Methods). Similarly, the orbital solution from radio-pulsar timing yields 
^psR = 30.008235 ± 0.000016 km s~^ for the pulsar. Combined, these constraints imply a mass 
ratio, q = Mpsr/Mwd = K^d/K^sr = H-VO + 0.13. 

Modeling of the Balmer-series lines in a high signal-to-noise average spectrum formed 
by the coherent addition of individual spectra (Fig. lb) shows that the WD has an effective 
temperature of T^g = (10120 + 47stat ± 90sys)K and a surface gravity of logiQ(^[cms"^]) = 
(6.035 + 0.032stat ± 0.060sys) dex. Here the systematic error is an overall estimate of uncer- 
tainties due to our fitting technique and flux calibration (8). We found no correlation of this 
measurement with orbital phase and no signs of rotationally-induced broadening in the spectral 
lines (8). Furthermore, we searched for variability using the ULTRACAM instrument (9) on the 
4.2-m William-Herschel Telescope at La Palma, Spain. The lightcurves, spanning 3 hours in to- 
tal, have a root- mean- square scatter of ~ 0.53, 0.07 and 0.08 mag in u', g' and / respectively 
and show no evidence for variability over the course of the observations. The phase-folded 
light-curve shows no variability either. Additionally, our calibrated magnitudes are consistent 
with the SDSS catalogue magnitudes, implying that the WD shone at a constant flux over this 
~ 5 yr timescale (8). 

Mass of the white dwarf The surface gravity of the WD scales with its mass and the inverse 
square of its radius (g = GMwd/^wd)- Thus, the observational constraints combined with a 
theoretical finite-temperature mass-radius relation for low-mass WDs yield a unique solution 
for the mass of the companion (10). Numerous such models exist in the literature, the most 
detailed of which are in good agreement for very low mass WDs (< 0.17 - 0.18 Mq), but differ 
substantially for higher masses [e.g. (11-13)]. The main reason for this is the difference in the 
predicted size of the hydrogen envelope, which determines whether the main energy source of 
the star is residual hydrogen burning (for "thick" envelopes) or the latent heat of the core (for 
"thin" envelopes). 

In the most widely accepted scenario, WDs lose their thick hydrogen envelope only if their 
mass exceeds a threshold. The exact location of the latter is still uncertain but estimated to 
be around ~ 0.17 - 0.22 M© [e.g. (11-13)]. Two other pulsars with WD companions, studied 
in the literature, strongly suggest that this transition threshold is indeed most likely close to 
0.2 Mq (10, 14). In particular, the WD companion of PSRJ1909-3744 has a well-determined 
mass of 0.20 Mq (15), a large characteristic age of a several Gyr and a WD companion that 
appears to be hot (10), suggesting that its envelope is thick. For this reason we base the WD 
mass estimate on cooling tracks with thick hydrogen atmospheres for masses up to 0.2 Mq, 
which we constructed using the "MESA" stellar evolution code (8, 16). Initial models were 
built for masses identical to the ones in (11) — for which previous comparisons have yielded 
good agreement with observations (14) — with the addition of tracks with 0.175 and 0.185 M© 



for finer coverage (Fig. 2). For masses up to 0.169Mo our models sliow excellent agreement 
with {11); our 0.196Mo model though is quite different, because it has a thick envelope instead 
of a think one. Being closer to the constraints for the WD companion to PSR J0348+0432, it 
yields a more conservative mass constraint: Mwd = 0.165 - 0.185 at 99.73% confidence (Fig. 3 
& Table 1), which we adopt. The corresponding radius is 7?wd = 0.046 - 0.092 R© at 99.73% 
confidence. Our models yield a cooling age of TcooI ~ 2 Gyr. 

Pulsar mass The derived WD mass and the observed mass ratio q imply a NS mass in the 
range 1.97-2.05 Mq at 68.27% or 1.90-2.18 M© at 99.73% confidence. Hence, PSRJ0348+0432 
is only the second NS with a precisely determined mass around 2Mo, after PSRJ1614-2230 
(2). It has a 3-cr lower mass limit 0.05 Mq higher than the latter, and therefore provides a verifi- 
cation, using a different method, of the constraints on the EOS of super-dense matter present in 
NS interiors (2, 17). For these masses and the known orbital period, GR predicts that the orbital 
period should decrease at the rate of P™ = (-2.58+°*") x 10"" s s"^ (68.27% confidence) due 
to energy loss through GW emission. 

Radio observations Since April 2011 we have been observing PSRJ0348-I-0432 with the 
1.4 GHz receiver of the 305-m radio telescope at the Arecibo Observatory, using its four Wide- 
band Pulsar Processors {18). In order to verify the Arecibo data, we have been independently 
timing PSRJ0348-I-0432 at 1.4 GHz using the 100-m radio telescope in Effelsberg, Germany. 
The two timing data sets produce consistent rotational models, providing added confidence 
in both. Combining the Arecibo and Effelsberg data with the initial GBT observations (7), 
we derive the timing solution presented in Table 1. To match the arrival times, the solution 
requires a significant measurement of orbital decay, P^ = (-2.73 + 0.45) x 10"" s s"^ (68.27% 
confidence). 

The total proper motion and distance estimate (Table 1) allows us to calculate the kinematic 
corrections to P\, from its motion in the Galaxy, plus any contribution due to possible variations 
of Newton's gravitational constant G: dP^, = 0.016 + 0.003 x 10" s s ^ This is negligible 
compared to the measurement uncertainty. Similarly, the small rate of rotational energy loss of 
the pulsar (Table 1) excludes any substantial contamination due to mass loss from the system; 
furthermore we can exclude substantial contributions to P^ from tidal effects (see {8) for details). 
Therefore, the observed P\, is caused by GW emission and its magnitude is entirely consistent 
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If we assume that GR is the correct theory of gravity, we can then derive the component 
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with the one predicted by GR: Pb/^b =1-05 + 0.18 (Fig. 3) 

m 
masses from the intersection of the regions allowed by q and P^ (Fig. 3): Mwd = 0.177^!!Q}g M, 



and MpsR = 2.07^[j2j Mq (68.27% confidence). These values are not too constraining yet. 
However, the uncertainty of the measurement of Pb decreases with T'^^^ (where T is the timing 
baseline); therefore this method will yield very precise mass measurements within a couple of 
years. 



Discussion 

PSR J0348+0432 as a testbed for gravity There are strong arguments for GR not to be valid 
beyond a (yet unknown) critical point, like its incompatibility with quantum theory and its 
prediction of the formation of spacetime singularities under generic conditions. Therefore, it 
remains an open question if GR is the final description of macroscopic gravity. This strongly 
motivates testing gravity regimes that have not been tested before, in particular regimes where 
gravity is strong and highly non-linear. Presently, binary pulsars provide the best high-precision 
experiments to probe strong-field deviations from GR and the best tests of the radiative prop- 
erties of gravity (19-23). Among these systems PSR J0348-I-0432 has a special role: it is the 
first massive (~ 2Mo) NS in a relativistic binary orbit. The orbital period of PSR JOB 4 8 -1-043 2 
is only 15 seconds longer than that of the double pulsar system, but it has ~ 2 times more 
fractional gravitational binding energy than each of the double pulsar NSs. This places it far 
outside the presently tested binding energy range [see Fig. 4a & (8)] . Because the magnitude 
of strong-field effects generally depends non-linearly on the binding energy, the measurement 
of orbital decay transforms the system into a gravitational laboratory for a previously untested 
regime, qualitatively very different from what was accessible in the past. 

In physically consistent and extensively studied alternatives, gravity is generally mediated 
by extra fields (e.g. scalar) in addition to the tensor field of GR (5). A dynamical coupling 
between matter and these extra fields can lead to prominent deviations from GR that only occur 
at the high gravitational binding energies of massive NSs. One of the prime examples is the 
strong-field scalarization discovered in (3). If GR is not valid, in the PSR J0348-I-0432 system 
where such an object is closely orbited by a weakly self-gravitating body, one generally ex- 
pects a violation of the strong equivalence principle that in turn leads to a modification in the 
emission of gravitational waves. While in GR the lowest source multipole that generates grav- 
itational radiation is the quadrupole, alternative gravity theories generally predict the presence 
of monopole and dipole radiation, on top of a modification of the other multipoles (5). For a 
binary system, the leading change in the orbital period is then given by the dipole contribution, 
which for a (nearly) circular orbit reads (8): 

Adipolar 4;r^G MpsrMwd . .2 /ix 

^b ^'T^M ri^ («PSR - aWD) , (1) 

C rb jWpsR -I- MwD 

where o-psR is the effective coupling strength between the NS and the ambient fields responsible 
for the dipole moment [e.g. scalar fields in scalar-tensor gravity], and q-wd is the same parameter 
for the WD companion. The WD companion to PSR J0348-I-0432 has a fractional gravitational 
binding energy (^grav/^wDf^) of just -1.2 x 10"^, and is therefore a weakly self-gravitating 
object. Consequently, q-wd is practically identical to the linear field-matter coupling ao, which 
is well constrained (laol < 0.004) in Solar System experiments (20, 24). 

For ttpsR, the situation is very different. Even if uq is vanishingly small, apsR can have 
values close to unity, due to a non-linear behavior of gravity in the interaction between matter 
and the gravitational fields in the strong-gravity regime inside NSs {3, 4). A significant apsR 



for NSs up to 1.47 M0 has been excluded by various binary pulsar experiments (8, 23). The 
consistency of the observed GW damping (Pb) with the GR predictions for PSR J0348+0432 
(Table 1) implies IccpsR - aol < 0.005 (95% confidence) and consequently excludes significant 
strong-field deviations, even for massive NSs of ~ 2 Mq. 

To demonstrate in some detail the implications of our results for possible strong-field devia- 
tions of gravity from Einstein's theory, we confront our limits on dipolar radiation with a specific 
class of scalar-tensor theories, in which gravity is mediated by a symmetric second-rank tensor 
field g* and by a long-range (massless) scalar field (p. Scalar-tensor theories are well motivated 
and consistent theories of gravity, extensively studied in the literature [e.g. {25, 26)]. For this 
reason, they are the most natural framework for us to illustrate the gravitational phenomena that 
can be probed with PSR J0348+0432. 

Concerning the EOS of NS matter, in our calculations we use the rather stifl" EOS ".20"of 
(27) that supports (in GR) NSs of up to 2.6 Mq. We make this choice for two reasons: i) a 
stifi'er EOS generally leads to more conservative limits when constraining alternative gravity 
theories, and ii) it is able to support even more massive NSs than PSR J0348-I-0432, which are 
likely to exist (28-30). Furthermore, in most of our conclusions a specific EOS is used only for 
illustrative purposes, and the obtained generic results are EOS independent. 

Fig. 4b illustrates how PSR J0348-I-0432 probes a non-linear regime of gravity that has not 
been tested before. A change in EOS and gravity theory would lead to a modified functional 
shape for the effective coupling strength, apsR. However, this would not change the general 
picture: even in the strong gravitational field of a 2 Mq NS gravity seems to be well described 
by GR and there is little space for any deviations, at least in the form of long-range fields, 
which influence the binary dynamics. Short range interactions, like massive Brans-Dicke grav- 
ity (31) with a sufficiently large scalar mass (heavier than ~ 10"^^eV/c^), cannot be excluded 
by PSR J0348 -1-0432. Nevertheless, as we will argue below, in combination with the upcom- 
ing ground-based GW detectors, this could lead to particularly illuminating insights into the 
properties of gravitational interaction. 

Constraints on the phase evolution of neutron-star mergers The first likely direct GW de- 
tection from astrophysical sources by ground-based laser interferometers, like the LIGO (Laser 
Interferometer Gravitational Wave Observatory, (32)) and the VIRGO (33) projects, will mark 
the beginning of a new era of GW astronomy (34). One of the most promising sources for these 
detectors are in-spiralling compact binaries, consisting of NSs and BHs, whose orbits are de- 
caying towards a final coalescence due to GW damping. While the signal sweeps in frequency 
through the detectors' typical sensitive bandwidth between about 20 Hz and a few kHz, the 
GW signal will be deeply buried in the broadband noise of the detectors (34). To detect it, one 
will have to apply a matched filtering technique, i.e. correlate the output of the detector with a 
template wave form. Consequently, it is crucial to know the binary's orbital phase with high 
accuracy for searching and analyzing the signals from in-spiraling compact binaries. Typically, 
one aims to lose less than one GW cycle in a signal with ~ 10"* cycles. For this reason, within 
GR such calculations have been conducted with great effort by various groups up to the 3.5 post- 



Newtonian order, i.e. all (non-vanishing) terms up to order (v/cy, providing sufficient accuracy 
for a detection (35-37). 

If the gravitational interaction between two compact masses is different from GR, the phase 
evolution over the last few thousand cycles, which fall into the bandwidth of the detectors, might 
be too different from the (GR) template in order to extract the signal from the noise. In scalar- 
tensor gravity for instance, the evolution of the phase is driven by radiation reaction, which is 
modified because the system loses energy to scalar GWs (38, 39). Depending on the difference 
between the effective scalar couplings of the two bodies, a a and ag, the 1.5 post-Newtonian 
dipolar contribution to the phase evolution could drive the GW signal many cycles away from 
the GR template. For this reason, it is desirable that potential deviations from GR in the inter- 
action of two compact objects can be tested and constrained prior to the start of the advanced 
GW detectors. For "canonical" 1.4 M© NSs and long-range gravitational fields, this has already 
been achieved to a high degree in binary pulsar experiments, e.g. {39). So far, the best con- 
straints on dipolar gravitational wave damping in compact binaries come from the observations 
of the millisecond pulsar PSR J1738-I-0333 {23). However, as discussed in detail above, these 
timing experiments are insensitive to strong-field deviations that might only become relevant in 
the strong gravitational fields associated with high-mass NSs. Consequently, the dynamics of 
a merger of a 2 Mq NS with a "canonical" NS or a BH might have a significant contribution 
from dipolar GWs. With our constraints on dipolar radiation damping from the timing observa- 
tions of PSRJ0348-I-0432, given above, we can already exclude a deviation of more than ~ 0.5 
cycles from the GR template during the observable in-spiral caused by additional long-range 
gravitational fields, for the whole range of NS masses observed in nature (see Fig. 5, and {8) 
for the details of the calculation). This compares to the precision of GR templates based on 
the 3.5 post- Newtonian approximation {35, 37). Furthermore, in an extension of the arguments 
in {38, 39) to massive NSs, our result implies that binary pulsar experiments are already more 
sensitive for testing such deviations than the upcoming advanced GW detectors. 

Finally, as mentioned before, our results on PSR J0348-I-0432 cannot exclude dipolar radi- 
ation from short-range fields. Hence, if the range of the additional field in the gravitational 
interaction happens to lie between the wavelength of the GWs of PSRJ0348-I-0432 and the 
wavelength of the merger signal (~ 10^ cm; ~ 10"^^ eV/c^), then the considerations concerning 
the applicability of the GR template given here do not apply. On the other hand, in such a case 
the combination of binary pulsar and LIGOA^IRGO experiments can be used to constrain the 
mass of this extra field. 

Formation, past and future evolution of the system The measured spin period P and spin- 
period derivative P of PSRJ0348-I-0432, combined with the masses and orbital period of the 
system (Table 1), form a peculiar set of parameters that gives insight to binary stellar evolu- 
tion. The short 2.46-hr orbital period is best understood from evolution via a common envelope 
where the NS is captured in the envelope of the WD progenitor, leading to efficient removal of 
orbital angular momentum on a short timescale of ~ 10^ yr {40). This implies that the NS was 
born with an initial mass close to its current mass of 2.01 M©, because very little accretion was 



possible. Whereas the slow spin period of ~ 39 ms and the unusually strong magnetic field (8) of 
a few 10^ G (Table 1) provide further support for this scenario, the low WD mass contradicts the 
standard common-envelope hypothesis by requiring a progenitor star mass smaller than 2.2 Mq, 
because more massive stars would leave behind more massive cores {8,41). For such low donor 
star masses, however, the mass ratio of the binary components is close to unity, leading to dy- 
namically stable mass transfer without forming a common envelope (42,43). One potential solu- 
tion to this mass discrepancy for common-envelope evolution is to assume that the original mass 
of the WD was > 0.4 Mq and that it was subsequently evaporated by the pulsar wind (44) when 
PSRJ0348 -1-0432 was young and energetic, right after its recycling phase (8). Such an evolution 
could also help explain the formation of another puzzling system, PSR J1744-3922 (45). How- 
ever, we find that this scenario is quite unlikely given that the observed spectrum of the WD in 
PSR J0348 -1-0432 only displays hydrogen lines, which is not expected if the WD was indeed a 
stripped remnant of a much more massive helium or carbon-oxygen WD. Furthermore, it is un- 
clear why this evaporation process should have come to a complete stop when the WD reached 
its current mass of 0.17 M©. A speculative hypothesis to circumvent the above-mentioned prob- 
lems would be a common-envelope evolution with hypercritical accretion, where ~ 0.6 M© of 
material was efficiently transferred to a 1.4 M© NS (8, 46). 

An alternative, and more promising, formation scenario is evolution via a close-orbit low- 
mass X-ray binary (LMXB) with a 1 .0 - 1.6 Mq donor star that suffered from loss of orbital 
angular momentum due to magnetic braking (43,47,48). This requires a finely tuned truncation 
of the mass-transfer process which is not yet understood in detail, but is also required for other 
known recycled pulsars (8) with short orbital periods of P\, < S hr and low-mass helium WD 
companions with Mwd ~ 0.14 - 0.18 Mq. The interplay between magnetic braking, angular 
momentum loss from stellar winds (possibly caused by irradiation) and mass ejected from the 
vicinity of the NS is poorly understood and current stellar evolution models have difficulties 
reproducing these binary pulsar systems. One issue is that the converging LMXBs most often 
do not detach but keep evolving with continuous mass transfer to more and more compact 
systems with P\, < Ihr and ultra- light donor masses smaller than 0.08 Mq. 

Using the Langer stellar evolution code (8), we have attempted to model the formation of 
the PSR J0348-I-0432 system via LMXB evolution (Fig. 6). To achieve this, we forced the donor 
star to detach its Roche lobe at Pb ~ 5 hr, such that the system subsequently shrinks in size to 
its present value of P^ =^ 2.46 hr due to GW radiation within 2 Gyr, the estimated cooling age 
of the WD. An illustration of the past and future evolution of PSR J0348-I-0432 from the two 
different formation channels is shown in Fig. 7. 

An abnormality of PSR J0348-I-0432 in view of the LMXB model is its slow spin period of 
P ~ 39 ms and, in particular, the high value for the spin period derivative, P = 2.41 x 10"^^ s s"^ 
These values correspond to an inferred surface magnetic flux density of 5 ~ 2 x 10^ G, which is 
high compared to most other recycled pulsars (49). However, a high B value naturally explains 
the slow spin period of PSR J0348-I-0432 from a combination of spin-down during the Roche- 
lobe decoupling phase (50) and subsequent magnetic dipole radiation from this high-magnetic- 
field pulsar (8, 49). Another intriguing question concerning this evolutionary channel is the 
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spread in NS masses. In the five currently known NS-WD systems with Pb < S hr, the NS 
masses span a large range of values, ranging from ~ 1.4 up to 2.0 Mq. The lower masses imply 
that the mass transfer during the LMXB phase is extremely inefficient — only about 30% of 
the material leaving the donor is accreted by the NS (14, 15). If this is indeed the case, and one 
assumes that the physical processes that lead to the formation of these systems are similar, it 
is likely that PSR J0348+0432 was born with an initial mass of 1.7 + 0.1 Mq, providing further 
support for a non-negligible fraction of NSs born massive (41). 

Emission of GWs will continue to shrink the orbit of PSRJ0348+0432 and in 400 Myr 
(when P\, =^ 23 min) the WD will fill its Roche lobe and possibly leave behind a planet orbiting 
the pulsar (51, 52). Alternatively, if PSR J0348+0432 is near the upper-mass limit for NSs then 
a BH might form via accretion-induced collapse of the massive NS in a cataclysmic, y-ray 
burst-like event (53). 

Materials & Methods 

Radial velocities and atmospheric parameters A detailed log of the VLT observations can 
be found in (8) (Figure SI & Table SI). We extracted the spectra following closely the method 
used in (14) and compared them with template spectra to measure the radial velocities. Our 
best fits for the WD had reduced ^^ minimum values of;^f^gjj^jj^ = 1.0-1.5 (8). Uncertainties 
were taken to be the difference in velocity over which ;^:'^ increases by;^^^^^^.^^ to account for the 
fact that;^f^gjj^jj^ is not equal to unity (14). After transforming the measurements to the reference 
frame of the Solar System Barycenter (SSB), we folded them using the radio-timing ephemeris 
described below. We then fitted for the semi- amplitude of the radial velocity modulation, ^wd^ 
and the systemic radial velocity with respect to the SSB, y, assuming a circular orbit and keeping 
the time of passage through the ascending node, Tasc, fixed to the best-fit value of the radio- 
timing ephemeris. Our solution yields ^wd = 351 + 4kms~' and y = -1 + 20kms~^ (8). 

Details of the Balmer lines in the average spectrum of PSR J0348 -1-0432, created by the co- 
herent addition of the individual spectra shifted to zero velocity, are shown in Fig. lb. We mod- 
eled the spectrum using a grid of detailed hydrogen atmospheres (54). These models incorpo- 
rate the improved treatment of pressure broadening of the absorption lines presented in (55). As 
mentioned above, our fit yields Tes = (10120 + 35stat ± 90sys) K for the effective temperature and 
logiQ^ = (6.042 + 0.032stat±0.060sys) for the surface gravity (5). The ;^f^ map shown in Fig. 2a is 
inflated to take into account systematic uncertainties. The average spectrum was also searched 
for rotational broadening. Using the analytic profile of (56) to convolve the model atmospheres, 
we scanned the grid of velocities < Vr sin / < 2000 km s"^ with a step size of 100 km s ^ The 
result is consistent with no rotation and our 1-cr upper limit is Vr sin / < 430 km s^ 

Modeling of the white-dwarf mass Low-mass WDs are thought to form naturally within 
the age of the Universe via mass transfer in a binary, either through Roche-lobe overflow or 
common-envelope evolution. In both cases, the WD forms when the envelope mass drops below 
a critical limit, which depends primarily on the mass of the stellar core, forcing the star to 



contract and detach from its Roche lobe. After the contraction, the mass of the relic envelope 
is fixed for a given core mass, but further reduction of its size may occur shortly before the 
star enters the final cooling branch due to hydrogen shell flashes which force the star to re- 
expand to giant dimensions. Additional mass removal via Roche-lobe overflow as well as rapid 
shell hydrogen burning through the CNO cycle may then lead to a decrease of the envelope 
size and aff'ect the cooling history and atmospheric parameters. To investigate the consequence 
of a reduced envelope size for the WD companion to PSRJ0348-I-0432, we constructed WD 
models in which we treat the envelope mass as a free parameter (8). For the WD companion 
to PSR J0348-I-0432, an envelope mass below the critical limit for hydrogen fusion is not likely 
for two main reasons: 

First, for a pure helium composition, the observed surface gravity translates to a WD mass of 
~ 0.15 M© and a cooling age of ~ 20 Myr, which is anomalously small. Such a small age would 
also imply a large increase in the birth and in-spiral rate of similar relativistic NS-WD systems 
(57). Furthermore, post-contraction flash episodes on the WD are not suflicient to remove the 
entire envelope. Therefore, creation of a pure helium WD requires large mass loss rates before 
the progenitor contracts, which is unlikely. For small progenitor masses (< 1.5 Mq) large mass 
loss prevents contraction and the star evolves to a semi-degenerate companion on a nuclear 
timescale that exceeds the age of the Universe. For more massive progenitors (> 1.5 Mq) the 
core grows beyond ~ 0.17 M© in a short timescale and ultimately leaves a too-massive WD. 

Second, even for envelope hydrogen fractions as low as Xavg = 10"^, the observed temper- 
ature and surface gravity cannot be explained simultaneously: The low surface gravity would 
again require a small mass of ~ 0.15 Mq. However, in this case the surface hydrogen acts 
like an insulator, preventing the heat of the core from reaching the stellar surface. As a result, 
temperatures as high as 10000 K can only be reached for masses above ~ 0.162 M©. 

Past a critical envelope mass, the pressure at the bottom of the envelope becomes high 
enough to initiate hydrogen-shell burning. The latter then becomes the dominant energy source 
and the evolutionary time-scale increases; the radius of the star grows by ~ 50% (depending on 
the mass), expanding further for larger envelopes. The dependence of the surface gravity on the 
radius implies that the observed value translates to a higher mass as the envelope mass increases. 
Therefore, the most conservative lower limit for the WD mass (and thus for PSR J0348 -1-0432, 
given the fixed mass ratio) is obtained if one considers models with the absolute minimum 
envelope mass required for hydrogen burning. In this scenario, the mass of the WD is in the 
range 0.162 - 0.181 M© at 99.73% confidence (8). Despite this constraint being marginally 
consistent with our observations (S), it is not likely correct due to the high degree of fine-tuning. 

For these reasons we have adopted the assumption that the WD companion to PSR J0348 -1-0432 
has a thick envelope as generally expected for WDs with such low surface gravity and high tem- 
perature. 

Radio-timing analysis The Arecibo observing setup (8) and data reduction are similar to the 
well tested ones described in (23). Special care is taken with saving raw search data, which al- 
lows for iterative improvement of the ephemeris and eliminates orbital-phase dependent smear- 
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ing of the pulse profiles, which might contaminate the measurement of P^ {58). From this 
analysis we derive 7773 independent measurements of pulse times of arrival (TOAs) with a 
root-mean-square (rms) uncertainty smaller than 10/is. Similarly the Effelsberg observations 
yield a total of 179 TOAs with uncertainties smaller than 20jus. 

We use the tempo2 timing package {59) to derive the timing solution presented in Table 1, 
using 8121 available TOAs from GBT (7), Arecibo and Effelsberg. The motion of the ra- 
diotelescopes relative to the barycenter of the Solar System was computed using the DE/LE 
421 Solar System ephemeris {60), published by the Jet Propulsion Laboratories. The orbit of 
PSRJ0348 -1-0432 has a very low eccentricity, therefore we use the "ELLl" orbital model {61) 
to describe the motion of the pulsar. 

For the best fit, the reduced x^ of the timing residuals (TOA minus model prediction) is 
1.66, a result similar to what is obtained in timing observations of other millisecond pulsars. 
The overall weighted residual rms is 4.6 /is. There are no unmodeled systematic trends in the 
residuals; either as a function of orbital phase or as a function of time. Therefore x^ > ^ 
is most likely produced by under-estimated TOA uncertainties. We increased our estimated 
TOA uncertainties for each telescope and receiver to produce a reduced x^ of unity on short 
timescales; for our dominant dataset (Arecibo) the errors were multiplied by a factor of 1.3. 

This produces more conservative estimates of the uncertainties of the timing parameters; 
these have been verified using the Monte Carlo statistical method described in {23): when all 
parameters are fitted, the Monte Carlo uncertainty ranges are very similar to those estimated by 
TEMPo2. As an example, tempo2 estimates P^ = (-2.73 ±0.45) x 10"^-' s s"' (68.27% confidence) 
and the Monte Carlo method yields P^ = (-2.72 + 0.45) x 10"'^ s s"^ (68.27% confidence), in 
excellent agreement. The observed orbital decay appears to be stable; no higher derivatives of 
the orbital period are detected {8). 
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Table 1. 

Observed and Derived Parameters for the PSR J0348+0432 system 

Timing parameters for the PSRJ0348+0432 system, indicated with their 1-cr uncertainties as 

derived by tempo2 where appropriate (numbers in parentheses refer to errors on the last digits). 

The timing parameters are calculated for the reference epoch MJD 56000, and are derived from 

TOAs in the range MJD 54872 - 56208. 

*For these timing parameters we have adopted the optically derived parameters (see text for 

details). 
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Optical Parameters 



Effective temperature, T^g (K) 10120 ± 47stat ± 90sys 

Surface gravity, log,o(^[cm s"']) 6.035 ± 0.032stat + 0.060sys 

Semi-amplitude of orbital radial velocity, /Two (km s"') 35 1 ± 4 

Systemic radial velocity relative to the Sun, y (km s"') . - 1 ± 20 

Timing Parameters 

Right ascension, a (J2000) 03" 48" 43^639000(4) 

Declination, 6 (J2000) +04° 32' 1 1'.'4580(2) 

Proper motion in right ascension, //„ (mas yr"') +4.04(16) 

Proper motion in declination, /ia (mas yr"^) +3.5(6) 

Parallax, n^ (mas) 0.47* 

Spin frequency, v (Hz) 25.5606361937675(4) 

First derivative of v, v (lO^^^ Hz s"') -0.15729(3) 

Dispersion measure, DM (cm^^ pc) 40.46313(1 1) 

First derivative of DM, DM1 (cm^^ pc yr^') -0.00069(14) 

Orbital period, Pt (d) 0.102424062722(7) 

Time of ascending node, T^sc (MJD) 56000.084771047(1 1) 

Projected semi-major axis of the pulsar orbit, X (It-s) .. 0.14097938(7) 

77 = esinw (+1.9 ± 1.0) x lO"*" 

K = ecosco (+1.4+ 1.0) X lO"*" 

First derivative of Pb, Pb (lO^'^ss"') -0.273(45) 

Derived Parameters 

Galactic longitude, / 183°3368 

Galactic latitude, b -36?7736 

Distance, d (kpc) 2. 1(2) 

Total proper motion, // (mas yr"') 5.3(4) 

Spin period, P (ms) 39.1226569017806(5) 

First derivative of P, P (10""^ s s"') 0.24073(4) 

Characteristic age, t^ (Gyr) 2.6 

Transverse magnetic field at the poles, Bq (10^ G) ~ 2 

Rate or rotational energy loss, E (10^^ erg s"') ~ 1.6 

Mass function, / (Mq) 0.000286778(4) 

Mass ratio, q = Mpsr/Mwd 1 1.70(13) 

White dwarf mass, Mwd (Mq) 0.172(3) 

Pulsar mass, MpsR (Mq) 2.01(4) 

"Range" parameter of Shapiro delay, r (/is) 0.84718* 

"Shape" parameter of Shapiro delay, s = sin / 0.64546* 

White dwarf radius, T^wd (Rg) 0.065(5) 

Orbital separation, a (10^ m) 0.832 

Orbital separation, a (Rg) 1 .20 

Orbital inchnation, / 40?2(6) 

Pb predicted by GR, P*^^ (lO^'^ s s"') -0.258+J}™^ 

Pb/Z'b'^ 1-05 ± 0.18 

Time until coalescence, t„, (Myr) ~ 400 
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Figure 1. 
Radial Velocities and Spectrum of the White Dwarf Companion to PSR J0348+0432. 
Upper: Radial velocities of the WD companion to PSR J0348+0432 plotted against the orbital 
phase (shown twice for clarity). Over-plotted is the best-fit orbit of the WD (blue line) and 
the mirror orbit of the pulsar (green). Lower: Details of the fit to the Balmer lines (Hfi to 
H12) in the average spectrum of the WD companion to PSRJ0348-I-0432 created by the co- 
herent addition of 26 individual spectra shifted to zero velocity. Lines from Hfi (bottom) to 
H12 are shown. The red solid lines are the best-fit atmospheric model (see text). Two models 
with (reff,logio^) = (9900 K, 5.70) and iT^fi■,\og^^^g) =(10200 K, 6.30), each ~ 3-ct off from 
the best-fit central value (including systematics) are shown for comparison (dashed blue lines). 
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Figure 2. 
Mass Measurement of the White Dwarf Companion to PSR J0348+0432 
Left: Constraints on effective temperature, Teff, and surface gravity, g, for the WD compan- 
ion to PSRJ0348+0432 compared with theoretical WD models. The shaded areas depict the 
X^ ~/fmin ~ •^••^' ^-2 ^'^^ 1^-^ intervals (equivalent to 1, 2 and 3-cr) of our fit to the aver- 
age spectrum. Dashed lines show the detailed theoretical cooling models of (ii). Contin- 
uous lines depict tracks with thick envelopes for masses up to ~ 0.2 Mg that yield the most 
conservative constraints for the mass of the WD. Right: Finite-temperature mass-radius re- 
lations for our models together with the constraints imposed from modeling of the spectrum 
(see text). Low mass - high temperature points are an extrapolation from lower temperatures. 
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Figure 3. 
System Masses and Orbital-Inclination Constraints 

Constraints on system masses and orbital inclination from radio and optical measurements of 
PSRJ0348+0432 and its WD companion. Each triplet of curves corresponds to the most likely 
value and standard deviations (68.27% confidence) of the respective parameters. Of these, two 
(the mass ratio q and the companion mass Mwd) are independent of specific gravity theories 
(in black). The contours contain the 68.27 and 95.45% of the two-dimensional probability 
distribution. The constraints from the measured intrinsic orbital decay (/'{,"', in orange) are 
calculated assuming that GR is the correct theory of gravity. All curves intersect in the same 
region, meaning that GR passes this radiative test (8). Left: cos /-Mwd plane. The gray region 
is excluded by the condition MpsR > 0. Right: Mpsr-Mwd plane. The gray region is excluded 
by the condition sin / < 1 . The lateral panels depict the one-dimensional probability-distribution 
function for the WD mass (right), pulsar mass (upper right) and inclination (upper left) based 
on the mass function, Mwd and q. 
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Figure 4. 
Probing Strong Field Gravity with PSR J0348+0432 

Left: Fractional gravitational binding energy as a function of the inertial mass of a NS in GR 
(blue curve). The dots indicate the NSs of relativistic NS-NS (in green) and NS-WD (in red) 
binary-pulsar systems currently used for precision gravity tests (8). Right: Effective scalar cou- 
pling as a function of the NS mass, in the "quadratic" scalar-tensor theory of (4). For the linear 
coupling of matter to the scalar field we have chosen ao = 10""*, a value well below the sensitiv- 
ity of any near- future Solar System experiment [e.g. GAIA (62)] . The solid curves correspond to 
stable NS configurations for different values of the quadratic coupling jSq: -5 to -4 (top to bot- 
tom) in steps of 0.1. The yellow area indicates the parameter space allowed by the best current 
limit on |apsR - orol (25), while only the green area is in agreement with the limit presented here. 
PSR J0348 -1-0432 probes deeper into the non-linear strong-field regime due to its high mass. 
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Figure 5. 
Constraints on the Phase Offset in Gravitational Wave Cycles in the LIGO/VIRGO bands 

Maximum offset in GW cycles in the LIGOA^IRGO band (20 Hz to a few kHz) between the GR 
template and the true phase evolution of the in- spiral in the presence of dipolar radiation, as a 
function of the effective coupling of the massive NS for two different system configurations: a 
2 Mq NS with a 1 .25 Mq NS (NS-NS), and a merger of a 2 Mq NS with a 10 Mq BH (NS-BH). 
In the NS-NS case, the green line is for as = ffo, and the gray dotted line represents the most 
conservative, rather unphysical, assumption oro = 0.004 and as = (8). In the NS-BH case, as 
is set to zero (from the assumption that the no-hair theorem holds). The blue line is for ao = 
0.004 (Solar System limit for scalar-tensor theories), and the purple line represents ao = 0. The 
gray area to the right of the red line is excluded by PSR J0348-I-0432. In this plot there is no 
assumption concerning the EOS. 
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Figure 6. 
Past and Future Orbital Evolution of PSR J0348+0432 

Formation of PSR J0348+0432 from our converging LMXB model calculation. The plot shows 
orbital period as a function of time (calibrated to present day). The progenitor detached from its 
Roche lobe about 2 Gyr ago (according to the estimated cooling age of the WD) when Pb - 5 hr, 
and since then GW damping reduced the orbital period to its present value of 2.46 hr (marked 
with a star). 
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Figure 7. 
Possible Formation Channels and Final Fate of PSR J0348+0432 

An illustration of the formation and evolution of PSR J0348+0432. The zero-age main sequence 
(ZAMS) mass of the NS progenitor is likely to be 20 - 25 Mq, whereas the WD progenitor had 
a mass of 1.0 - 1.6 Mq (LMXB) or 2.2 - SM© (common envelope, CE), depending on its 
formation channel. In ~ 400 Myr (when Pb - 23 min) the WD will fill its Roche lobe and the 
system becomes an ultra-compact X-ray binary (UCXB) leading to the formation of a BH or a 
pulsar with a planet. 
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Supporting Online Material 

VLT spectral observations and analysis We observed the companion of PSR J0348+0432 
during December 19 & 20, 201 1 with the F0RS2 (63) instrument on Unite-Telescope 1 (Antu), 
using its blue sensitive E2V CCD detectors and the G1200B grism. This setup delivers a res- 
olution of 0.36 Aperbinned-by-two-pixel along dispersion and O'.'25perbinned-by-two-pixel 
along the spatial direction. Because of the short orbital period of the binary, we chose a rela- 
tively wide 1" slit to avoid severe radial velocity smearing (by reducing the exposure time) and 
minimize possible dispersion losses not corrected by the dispersion corrector of the instrument. 
However, this choice may potentially result in systematic offsets in radial velocity measure- 
ments due to non-uniform illumination of the slit. To monitor these effects we rotated the slit 
by 134° 8 (north-through east) with respect to the parallactic angle to include a bright nearby 
star for local flux and velocity calibration (Fig. SI). Our setup covers the spectral range from 
~ 3700 to 5200 A with a resolution ranging from ~ 2 to 3 A depending on the seeing. During 
the first night the conditions were good to photometric and the seeing varied between ~ (X.'V 
and r.'2. The second night was sporadically plagued with thin cirrus and the seeing ranged 
from ~ (y.'9 to I'.'V. Bias, flat and Mercury-Cadmium (HgCd) frames for wavelength calibration 
were collected during day-time after each run. We collected a total of 34 spectra of the white 
dwarf companion to PSR J0348 -1-0432 and the nearby comparison star. Of these, 22 had 800-s 
exposures and were taken with the slit rotated by the angle mentioned above, 4 were taken with 
850-s exposures during bad weather instances and 8 with the slit rotated by a slightly different 
angle during experimental stages. In addition, we collected 2 spectra of the comparison through 
a wider, 2'.' 5 slit and spectra of several flux standards at the beginning and the end of each run 
through both 1" and 2'.'5 slits. 

We reduced the data using routines inside the Munich Image and Data Analysis System 
(midas). Our analysis, from cosmetic corrections to extraction of spectra, is identical to that 
followed for the white dwarf companion to PSRJ1738-I-0333 and is described in detail else- 
where (14). The dispersion solution has root-mean- square (rms) residuals of ~ 0.03 A for 18 
lines. Flux calibration was performed separately for each night by comparing flux-standard ob- 
servations (Table SI) with high S/N templates or appropriate white dwarf model spectra (54). 
Overall, the response curves from each standard are consistent with each other, with the largest 
differences (up to 10%) observed at short wavelengths (A < 4000 A); we use their average for 
flux calibration. Prior to the latter, we corrected the narrow-slit spectra for wavelength depen- 
dent slit losses using the wide slit spectra of the comparison and accounted for atmospheric 
extinction using the average extinction curve for La Silla. 

Radial velocities We extracted the radial velocities of the white dwarf companion and the 
nearby comparison star following the procedure described in (14). First, we identified the 
nearby comparison star as being a type GIV star (with an uncertainty of about 2 subtypes) 
and used a high-resolution spectrum of the similar star HD 20807 (64) as a template. For the 
PSR J0348 -1-0432 white-dwarf companion we fitted a high S/N spectrum with a grid of DA 
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model atmospheres (54), used the best-fit template to measure radial velocities, averaged the 
zero-velocity spectra and finally re-fitted the average spectrum to determine the final template. 
We scanned a grid of velocities from -800 to -1-800 km s"^ with a step-size of 5 km s ^ The best 
fits had ;^:'^gj j^jj^ = 1-1.5 andx^^^^j^^j^ = 1-1 - 3.0 for the comparison star. As described in the 
main paper, we scaled the errors to account for the fact that^^^^j^j^^ was not equal to unity. 

The velocities of the comparison star show a peak-to-peak variation of ~ 40kms^\ much 
higher than the typical 0.8 km s"^ measurement error. While we find no evidence for binarity, 
the measurements form 5 distinct groups, each of which display a variability only marginally 
higher than the formal errors. These coincide with blocks of observations interrupted for target 
repositioning. The scatter of velocities is therefore clearly related to the instrument and most 
probably associated with positioning uncertainties. For this reason we chose to use velocities 
relative to the comparison. 

The best-fit solution using all available (barycentred) data gave ^wd = 346 + 6kms~^ and 
a systemic velocity of 6y = -i-8 + 4km s"' relative to the comparison with ;^^^gjj ^^.^^ = 2.78 for 
32 degrees of freedom (dof). However, 8 of the observations used here were taken with the 
slit at a different angle and the white dwarf's velocity relative to the comparison is thus most 
likely contaminated with an extra systematic shift due to slit rotation. For this reason we neglect 
these data. Using the homogeneous set of observations only and further rejecting one outlier 
with spuriously shaped continuum (no. 15 in Table SI) we obtain ^wd = 345 ± 4km s"^ and 
6y = +23 + 5kms"^ respectively with^^^^j^j^^ = 0.99 for 23 degrees of freedom (Fig. la). 

After correcting for the small effect of orbital smearing (sin(7r(?exp)/^b)/(^(4xp)/^b) = 
0.98636) we find a semi- amplitude of ^wd = 351 ± 4km s^^ The best-fit systemic velocity 
of PSRJ0348-I-0432 using the raw white dwarf velocity measurements is y = -1 + 6kms~^ 
Given the large scatter of the comparison's velocity we adopt y = -1 + 20km s"' with the 
uncertainty being a conservative estimate based on the scatter of the data. 

Average spectrum and atmospheric parameters We scanned a grid of models covering 
effective temperatures from Teff = 8000 to 25000 K with a step-size of 250 K and surface 
gravities ranging from log^ = 5.00 to log^ = 8.00 with a step-size of 0.25 dex. At each 
point of the grid we fitted for the normalization using a polynomial function of the wave- 
length to account for non-perfect flux calibration. Analysis of statistical errors is again identical 
to that followed for the white dwarf companion of PSRJ1738-I-0333 {14). We achieved the 
best fit to higher Balmer lines when excluding the continuum regions between 4000 - 4050, 
4180 - 4270 and 4400 - 4790 A, which had small irregularities due to leftover detector imper- 
fections: Teff = 10120+35 K and log g = 6.042+0.032 (Icr) with ;t^2^^^.^ = 1.02. To estimate the 
influence of systematics we varied the degree of the polynomial used for normalization (1st to 
5th degree), the spectral regions used for the fit (lines-only to whole spectrum) and the assumed 
spectral resolution (by steps of ~ 5%). We also searched for velocity smearing by checking 
the consistency of the solution in an average of spectra taken close to orbital conjunction and 
an average of spectra taken close to the nodes. Finally, we fitted each line (from HyS to H12) 
separately to verify the consistency of the fit over the spectrum and examined the influence 
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of our flux calibration by fitting the average uncalibrated spectrum. Overall, all tests gave fits 
consistent within statistical errors with only few exceptions that had (higher) central values that 
diff"ered by 120 K and 0.1 1 dex compared to the numbers above. The good agreement is prob- 
ably due to the high S/N of the spectrum. The values adopted in the main paper are based on 
the solution using a third degree polynomial and the systematic error is a conservative estimate 
based on the scatter of the different fits mentioned above. 

Spectroscopic modeling and the "high log g" problem Spectroscopic modeling of the Balmer 
lines in higher mass white dwarfs shows a spurious increase in surface gravity for stars with 
temperatures between ~ 8000 and 1 1000 K. This well-known problem is linked to the incom- 
plete treatment of convection in 1-D atmospheric models and disappears with the use of 3-D 
model atmospheres {65, 66). However, our modeling below shows that for the parameter space 
relevant to the PS R JOB 4 8 -1-043 2 companion, the atmosphere is not yet convective (e.g. convec- 
tion sets in at T^ff < 9300 K for Mwd = 0.17 Mq). Therefore this problem is very unlikely to be 
relevant for the mass determination presented here. 

Initial white-dwarf models To construct the white dwarf models presented in the main pa- 
per, we evolved solar composition stars (metal mass-fraction of Z = 0.02) with masses between 
1 .0 and 1 .5 Mq and applied a large mass-loss wind at various points on the Red Giant Branch 
(RGB). To constrain the upper limit of the envelope mass expected from natural binary evolu- 
tion, we removed the mass before the star enters the asymptotic RGB, letting the star evolve 
and contract naturally to become a white dwarf. Our upper limits agree well with the results 
of previous studies (11, 12, 67). Finally, to fully control the envelope mass of the white dwarf 
at the final stages of evolution we neglected hydrogen fusion through the CNO bi-cycle that is 
responsible for the hydrogen shell flashe^ 

In Fig. S2 we show the post-contraction white dwarf cooling age when T^s = 10000 K, as 
a function of the total hydrogen mass (after cessation of the mass transfer), for masses ranging 
from 0.155 to 0.185 Mq. For low envelope masses, hydrogen burning cannot be initiated and 
the white dwarf quickly radiates the latent thermal energy of the core and cools in a few Myr. 
The thick-envelope modes presented in the main text were constructed as above. 

Metallicity The metallicity of the white dwarf plays an important role in both regulating the 
CNO luminosity and changing the chemical profile of the stellar envelope. Qualitatively, our 
main models described above are in good agreement with the Z = 0.001 models of {69) for 
the parameter space relevant to the white dwarf companion to PSRJ0348 -1-0432. Specifically, 
their 0.172Mo track has a thick envelope and predicts a surface gravity of log^ = 6.13 for 
Teff = lOOOOK which is reached at a cooling age of Tcooi = 2.85 Gyr. This agreement is not 
surprising given that CNO burning is neglected in our analysis, convective mixing has not yet set 



'For a white dwarf at the final cooUng branch, CNO luminosity accounts for less than 5% of the total energy 
budget. Hence, it is safe to neglect it without influencing the macroscopic characteristics of the models {68). 
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in at r = 10000 K and consequently metals are absent from surface layers due to gravitational 
settling. Therefore we consider that any uncertainties due to metallicity are small and anyway 
included in our adopted errors. 

Input physics of the stellar evolution models Stellar models used in our analysis were con- 
structed using the 1-D stellar evolution code "star" provided with the Modules for Experi- 
ments in Stellar Astrophysics (mesa) (16). star solves for the equations of hydrostatic equi- 
librium, nuclear energy generation, convection and time-dependent element diffusion using a 
self-adaptive non-Lagrangian mesh and analytic Jacobians. We used default options for the 
equation-of-state, radiative and neutrino opacities, thermonuclear and weak reaction rates de- 
scribed in (16) and references therein. We implemented the mixing length theory of convection 
from (70) that takes into account radiative losses near the outer layers of the star. Diffusion 
was taken into account using the method and coefficients from (71) and transport of material 
was calculated using the method described in (72) after grouping the elements in "classes" in 
terms of atomic mass ranges. Finally, boundary atmospheric conditions were calculated using 
the gray-atmosphere approach of (73). 

Photometry A photometric campaign on the white dwarf companion to PS R JOB 4 8 -1-043 2 
was carried out during February 1, 2012 using the ULTRACAM instrument (9) on the 4.2-m 
William-Herschel Telescope at La Palma, Spain. The data were reduced using the standard 
ULTRACAM pipeline (Fig. S3). 

The lightcurves have an rms scatter of ~0.53, 0.07 and 0.08 mag in u', g' and r' respectively 
and show no evidence for variability over the course of the observations. The phase-folded 
light-curve shows no variability either. Additionally, our calibrated magnitudes are consistent 
with the SDSS catalogue magnitudes implying no significant variability at the ~ 5 yr time-scale. 

Qualitatively, this result supports the use of the PSRJ0348-I-0432 system as a gravitational 
laboratory (see timing analysis below). In what follows, we discuss the limits on various param- 
eters in more detail. Three effects that cause phase-dependent variability are: deformation of 
the white dwarf by tides raised by the neutron star, irradiation by the pulsar wind, and Doppler 
boosting caused by the white dwarf's orbital motion. For a circular orbit, the combined modu- 
lation in photon rate n^ is given by 

Ariy/nyfi ^ /ell ( -^ J q sin^ i cos(47r0) + /db — ^ sin / cos(2n(/)) - /in-TTT^ sin / sinilncp), (2) 

with /ell, /db, and /„ factors of order unity describing the observability in a given filter, < < 1 
the orbital phase, and T^^ = Lp^JA-na^cr =^ 2400 K the effective temperature corresponding to 
the pulsar flux incident on the white dwarf. We find that all terms should be small. For the tidal 
deformation, /en = -3(15 -I- mi)(1 -I- ti)/20(3 - ui) = 1.75, where we use linear approximations 
for limb and gravity darkening, with coefficients ui = 0.36 (74) and ti = 1 (appropriate for a 
radiative atmosphere). Thus, the expected modulation is 1.5 x 10^. For the Doppler boosting, 
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approximating the white dwarf as a black-body emitter, fn, - a exp Qr/(exp a - I) - 2.6, where 
a = hc/AkT^ff =^ 2.8 (75), with A - 550 nm the typical observing wavelength. Hence, the 
expected amplitude is ~ 3 x 10"-'. Finally, for the irradiation, /irr = (1 - A)/db ^ 2, where the 
maximum is for albedo A ^ 0. Thus, irradiation could cause a modulation of up to ~ 1 .2 x 10 '*. 

Fitting the observed lightcurves with a function of the form Any/{ny) = 1 + aeii cos(4;7r0) + 
adb cos(2:;r0) - (3,>r sin(2:;r0), we find good fits (x^^^ - 1) but no significant detections, with 
averaged amplitudes of the higher S/N r and g band lightcurves of agii = 0.003 + 0.003, a^, = 
0.003 ± 0.003 and a^^ = 0.006 ± 0.004. The marginal irradiation signal would correspond to 
a temperature difference between the irradiated and non-irradiated side of ~ 100 K, which is 
substantially larger than the expected difference of 2 K. Even if confirmed, however, this would 
not affect our inferred radial velocity amplitude or white dwarf parameters. 

Finally, another possible source of variability is quadrupole moment variations of the white 
dwarf (68): these typically change the star's luminosity by a few per-cent (e.g. ~ 20% for the 
only three known cases of pulsating low-mass white dwarfs (76, 77)) and result in changes of 
the orbital period P^ through classical spin-orbit coupling (78). To our knowledge, all possible 
mechanisms for such variations would result in modulations much higher than the precision 
of our lightcurve. Therefore we can neglect this effect and assume that the star is in equilib- 
rium. Our assumption is further supported by the lack of second or higher-order derivatives 
in the measured orbital period (see below) and recent theoretical findings (79) that locate the 
instability strip for g-mode oscillations outside the parameter space relevant for the white dwarf 
companion to PSR J0348+0432. 

The SDSS photometry places a constraint on the distance to the system. Adopting the 
model of (80) for the interstellar reddening and the 0. 169 Mq cooling track of (11), we find that 
the luminosities (Fig. S3) are consistent with a distance of d - 2.1 kpc (and a reddening of 
Ay ~ 0.7). Given the uncertainties in the models the error is difficult to estimate but it should be 
better than ~ 10%. Our estimate is also consistent with the distance of Jdm ~ 2 kpc implied by 
the dispersion measure (DM) of the pulsar and the NE2001 model for the Galactic free electron 
density (81). 

Radio observations The observing setup for the Arecibo telescope is identical to the well- 
tested setup described in (23), with the exception of one WAPP now being centered at 1610 MHz 
instead of 1310 MHz; the former band is cleaner and its use improves the precision of our DM 
measurements. Also, as in the former case, data are taken in search mode and processed off-line. 
This allows for iterative improvement of the pulsar ephemeris which is important at the early 
stages when the timing parameters are not yet very precise. With each improved ephemeris, 
we de-disperse and re-fold the data, obtaining pulse profiles with higher S/N that yield more 
accurate pulse times-of-arrival (TOAs). This helps to avoid orbital-phase dependent smearing 
and timing artefacts, which may corrupt the determination of orbital parameters, particularly 
the orbital phase and orbital period variation (58). 

We de-disperse and fold the radio spectra following the procedure described in (23). TOAs 
are derived every 4 minutes to preserve the orbital information in the signal. The pulse profile 
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template, resulting from more than 1 hour of data, is displayed in Fig. |4j Although the pulse 
profile changes significantly from 350 to 2200 MHz (7), the changes within the band of the L- 
wide receiver used for timing (1 100-1660 MHz, also displayed in Fig. |4]) are small enough for 
us to consider this single average profile taken at 1410 MHz as a good template for all the data. 
The latter is cross-correlated with every 4-minute/25 MHz- wide pulse profile in the Fourier 
domain (82, 83) and the phase offset that yields the best match is used to derive the topocentric 
TOA of a reference sub-pulse (normally that closest to the start of each sub-integration). The 
results described below are obtained using 7773 TOAs with stated rms uncertainty smaller than 
lOjus. 

In order to verify the Arecibo data we have been timing PSRJ0348 -1-0432 with the 100-m 
radio telescope in Effelsberg, Germany, which has a very different observing system. The po- 
larization characterization of the radio emission of PSR J0348-I-0432, displayed in the top plot 
of Fig.|4| was made with this telescope. Overlaid on the polarization data is a theoretical Rotat- 
ing Vector Model (RVM). It is generally difficult to fit a RVM model to polarization data from 
recycled pulsars, but for PSR J0348 -1-0432 this model works surprisingly well. For instance, 
as explained in (82), the covariance between the angle between the spin and magnetic axis, a, 
and the angle between the spin axis and the line of sight ^, allows for a wide range of possible 
solutions (Fig. S5). However, if we assume that during the accretion episode that recycled the 
pulsar the spin axis of the pulsar was aligned with the orbital angular momentum (which has 
an angle / = 40:2 + 0°6 to the line of sight) then a - 45°. The minimum angle between the 
magnetic axis and the line of sight is then given hy /3 = ^ - a = -5°. 

Apart from the polarimetry, the Effelsberg data yielded a total of 179 high-quality TOAs. As 
can be seen in Fig. [6} these follow the Arecibo timing very closely, providing added confidence 
in both. 

Timing analysis The combined timing dataset contains 8121 TOAs. The TOA residuals ob- 
tained with the best ephemeris (Table 1) are displayed as a function of time in the top panel and 
as a function of orbital phase in the bottom plot of Fig. [6} To derive the ephemeris in Table 1 
(using TEMPo2) we increased the TOA uncertainties by factors of 1.3 for the GBT and Arecibo 
data and by 1.8 for the Effelsberg data. This results in the residuals of each dataset having a 
normalized ;^^ of 1. Using these slightly increased (but more realistic) TOA uncertainties results 
in more conservative (i.e. larger) uncertainties for the fitted timing parameters. Globally, the 
residualshavea weighted rms of 4.6 /is and the reduced ;^f^ is 1.019 for 8102 degrees of freedom. 
The TOA uncertainties presented in Fig.|6]are those used to derive the timing solution. 

The orbit of PSR J0348 -1-0432 has a very low eccentricity, therefore we use the "ELLl" 
orbital model (61) to parametrize it jj This parametrization yields Keplerian and post-Keplerian 
parameters very weakly correlated with each other. In order to estimate the intrinsic ("real") 



^The ELLl timing model as implemented in the tempo2 software package is a modification of the DD timing 
model (84, 85) adapted to low-eccentricity binary pulsars. In terms of post-Keplerian observables, it contains all 
those which are numerically relevant for systems with e <K 1 . The "Einstein delay" term is not relevant for such 
systems and is therefore not taken into account. 
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eccentricity of the binary (Table 1) we adopt Mwd = 0.172 Mq and / = 4(K2 obtained from 
the optical observations. This assumption is safe because GR is known to provide a sufficiently 
accurate description of spacetime around weakly self-gravitating objects (86). According to 
(87), the orthometric amplitude of the Shapiro delay (which quantifies the time amplitude of 
the measurable part of the Shapiro delay) is hi, = 42 ns. Fitting for this quantity we obtain 
/?3 = 69 + 53 ns. This is 1-cr consistent with the prediction but the low relative precision of this 
measurement implies that we cannot determine Mwd and sin / independently from the existing 
timing data. A precise measurement of the component masses of this system from Shapiro delay 
would require an improvement in timing precision that is much beyond our current capabilities. 

Intrinsic orbital decay As described in the main text, we detect an orbital decay consistent 
with the prediction of General Relativity. When we say that this decay is stable, we mean that 
we detect no higher-order variations of the orbital frequency /b = 1 jPh nor large variations in 
X = ap sin i/c: 



d'f, 



^ -- -4.5 + 4.4xlO"2^Hzs"^ (3) 



= +7.4 + 4.4xl0"^^ss"\ (5) 



—J = +4.1 + 2.5xlO-^^Hzs-^ (4) 

dx 
dt 

where the values and 1-cr uncertainties were obtained using the tempo implementation of the 
BTX orbital model. In systems where the quadrupole moment of the white dwarf changes, we 
should expect such timing effects (88) plus significant photometric variations with orbital phase 
(discussed above). Since none are observed, the companion to PSRJ0348 -1-0432 is very likely 
to have a stable quadrupole moment. 

The constraints on the total proper motion // combined with the optically derived distance 
d = 2.1 + 0.2 kpc allow us to calculate the two kinematic corrections to the observed P\,. The 
more important one is the Shklovskii effect (89): 

pshk ^ p^l^ ^ 0.0U911ZI X 10"'' s s-\ (6) 

where we have adopted the 10% error-estimate on the distance. The second correction is caused 
by the difference of Galactic accelerations between the binary and the Solar System. Using the 
detailed procedure outlined in (23), we obtain: 

pAcc ^ p^^ ^ O.OOSrJ^Z' X 10"^' s s-i. (7) 

A third correction could arise from a possible variation of the gravitational constant G. Conser- 
vative limits are given by (90-93): 

P^ = -2P^% = (0.0003 + 0.0018) X 10" s s^ (8) 

G 
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where we used the latest limit on G from Lunar Laser Ranging (94). 

Adding these corrections, we obtain a total of ~ (+1.6 + 0.3) x 10"^^ s s"\ or about 0.006 
of the measured value. This is much smaller than the current measurement uncertainty and 
therefore we can conclude that, at the current precision limit, the observed value is intrinsic to 
the system. Its magnitude is entirely consistent with the GR prediction for the orbital decay 
caused by emission of gravitational waves: P\,/P^^ = 1.05 + 0.18. This agreement is depicted 
graphically in a cos / - Mwd and MpsR - Mwd diagram (Fig. 3); the consequences are discussed 
in the main paper and in detail further below. 

Mass loss contribution to P^ If the system is losing mass, that should cause a change in the 
orbital period (97): 

P* = 2'^K. (9) 

where Mj = MpsR + ^wd is the change of mass of both components. 

We now estimate both mass loss terms. The pulsar is losing rotational energy at a rate given 
by £■ = AnlpsRpP'^ = 1.6 x lO"*^ erg s~\ where /psR is the pulsar's moment of inertia, normally 
assumed to be lO'^^gcm^. This dominates the mass loss for the pulsar (91): 

—^ = =4.1xl0'2^s~\ (10) 

Most of this energy is emitted as a wind of relativistic particles, which we assume to be isotropic 
to first order. A fraction of this energy F = R^^/4a^ = 0.00074 (where a = xc(q +1)1 sin i = 
8.32 X 10^ m is the separation between components) strikes the surface of the white dwarf. 
This is the energy available to power mass loss from the white dwarf. Conservation of energy 
requires that 

EF = ^M^^v\ (11) 

where v is the velocity of the escaping particles. This equation shows that M increases as v 
decreases, however v must be at least equal to the escape velocity for the star to lose mass, i.e., 
v^/2 > GMwd/7?wd- Putting all the constraints together, we obtain: 

^ < 5.4 X 10-21 s-^ (12) 

Therefore, Mj =^ Mwd- Evaluating eqM we obtain Pf < 0.4 x 10"^^, which is ~ 5 x 10^ times 
smaller than the current uncertainty in the measurement of Pb. 

Tidal contribution to Pb We now calculate the orbital decay caused by tides. If these change 
the angular velocity of the white dwarf Qwd^ this will be compensated by a change in the orbital 
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period of the system P^. We can relate the two because of conservation of angular momentum: 

• T _ 3fcQwD / ^wD^bSin/ \^ 1 
'' 2nq(q + 1) \ xc ) r/ 

where r^ = -Qwd/^wd is the synchronization timescale and k = /wd/(^wd^wd)' where /wd 
is the white dwarf moment of inertia. For idealized white dwarfs (particularly those with a 
mass much below the Chandrasekhar limit) sustained solely by degeneracy pressure of non- 
relativistic electrons, apolytropic sphere with n = 1.5 provides a good approximation. For such 
stars, we have k = 0.2 (95). However, for this light white dwarf only the core is degenerate, and 
is surrounded by a deep non-degenerate layer that accounts for only about 5% of the mass of 
the star. Therefore, the mass distribution is much more centrally condensed than for ann = 1.5 
polytrope and the moment of inertia is much smaller. We therefore use the output of our white 
dwarf model calculations (see Fig 2. and above) to estimate that factor. For the model closer 
to the mean of the white dwarf mass distribution, with Mwd = 0.169Mo, ^wd = 0.069 R© and 
Tes = 9950 K we obtain k = 0.0267. We adopt this value in subsequent calculations. 

The only unknown parameters in this expression are Qwd and r.,. If r^ were much smaller 
than the characteristic age of the pulsar Tc = 2.6 Gyr (which is similar to the cooling age of the 
white dwarf, i.e., this number is likely to be a good approximation to the true age of the system), 
then the white dwarf rotation would already be synchronized with the orbit {Q.wd = '^.n/Pb). In 
this case the orbital decay would be slightly affected because, as the orbital period decreases, 
the white dwarf spin period would decrease at exactly the same rate in order to preserve tidal 
locking. The resulting exchange of angular momentum would change the orbital decay by a 
factor APb given by the ratio of the moment of inertia of the white dwarf and the binary: 

AP, /wD k /^wDsinA^^^^xlO-^ (14) 



PGR 4 ^(^+1)\ xc 

This means that, were the system synchronized, APb would be an insignificant correction given 
our current measurement precision. 

If the white dwarf is not yet synchronized, then r^ > Tc- In this case IIwd can be much 
larger than In/P^, but it must still be smaller than the break-up angular velocity, i.e., Qwd < 
(GMwd/^wd)'''^ = 0.0142 rads"'. These conditions for Qwd and t, yield Pj^ < 4.2x 10"^^ ss~^ 
Thus, even if the white dwarf were rotating near break-up velocity, P^ would still be two orders 
of magnitude smaller than the uncertainty in the measurement of P\,. We note, however, that the 
progenitor of the white dwarf was very likely synchronized with the orbit at formation, which 
had a period of ~ 5 hours (see below). When the white dwarf formed, fall-back of material 
within the Roche lobe into it would have spun it up, but not by more than 1 order of magnitude 
(e.g.. Appendix B2. 2 of (9(5)). Therefore, at formation Qwd wasof the order of 3.5x10"^ rads"'; 
this would yield Pj^ < 1.0 x 10"^^ ss^ 

Constraints on dipolar radiation and Scalar-Tensor gravity In scalar-tensor gravity, like 
for most other alternatives to GR, the dominant contribution to the GW damping of the orbital 
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motion of a binary system would come from the scalar dipolar waves, proportional to (oa - 
asf, where Qa and Qb denote the effective scalar-coupling constants of the two masses m^ 
and Mb, respectively, of the binary system. Such deviations should then become apparent as 
a modification in the orbital period decay observed in binary pulsars. In GR the emission of 
quadrupolar tensor waves enters the orbital dynamics at the 2.5 post-Newtonian (pN) level, 
which corresponds to corrections of order (v/c)^ in the equations of motion, v being a typical 
orbital velocity. A contribution from dipolar GWs enters already at the 1.5pN level, i.e. terms 
of order (v/c)^. As an example, in scalar-tensor gravity the change in angular orbital frequency 
Rb = 27r//'b for a circular orbit caused by gravitational wave damping up to 2.5pN order is given 
hy {39,97) 



nl 1+ QaUb 



96 /v\5 

K +(aA- 
5 \cl 


\cl 



(15) 



where 

V = [Gh.(1 -I- aAaBXrtiA + mB)nb]'^^ , (16) 

with Gh. denoting the bare gravitational constant, and Xa = mAliniA + m^) and Xb = niB/iniA + 
nis). The quantity /c, where a- = 1 in GR, holds terms arising from the emission of scalar 
quadrupolar waves and higher order terms of the scalar dipolar emission (39): 

K = I + - (qaXb + obXa)^ + di (ua - Ob) + ^2 (ctA - ^s)^ (IV) 



an 
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GR is recovered for G^. = G and qa = aB = 0. Equation ( |T5) ) can directly be confronted 
with the results compiled in Table 1, in combination with the p.d.f. of the white dwarf mass 
in Fig. 3, where we use index A for the pulsar and index B for the white dwarf companion|j 
One finds from the mass ratio q that Xa = q/(q + 1) = 0.9213 + 0.0008 and Xb = \/(q + \) = 
0.0787 + 0.0008. Furthermore, since G,{mA + mb) - GniBiq + 1), one has v/c = (0.001970 + 
0.000016) X (1 -I- qtaQ'b)'^^. With the observed change in the orbital frequency h^, = -InPb/Pl ~ 
(2.23 + 0.36) X 10~^°, which agrees with GR, one can infer the following constraint on (a a - ao) 



using equation (15) 



\aA -ao\< 0.005 (95% C.L.) . (20) 



^To our knowledge, the quantity d2 has been calculated here for the first time. 

^Strictly speaking, when using the masses of Table 1 in equation ( 15 1 one has to keep in mind the difference 



between the bare gravitational constant G* and Newton's gravitational constant G - G,(l + cq)^ as measured in a 
Cavendish-type experiment. However, since a^ < 10"^, we can ignore this difference in our calculations. 
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Our detailed calculations show that this limit is solely enforced by the dominant 1.5pN term of 
equation ( [T5] ) (see also equation (1) in the main text), and is practically insensitive to the values 
assumed by /3a and /3b- Consequently, as in {23, 98) the limit ( [20] ) can be seen as a generic limit, 
that is independent of the EOS. 

To illustrate how PSR J0348+0432 probes a new gravity regime, we present detailed calcu- 
lations based on a specific EOS and a specific class of alternative gravity theories. As an EOS 
we use the rather stiff EOS ".20" of {27), which supports (in GR) neutron stars of up to 2.6 Mq. 
Concerning the alternative gravity theories, we use the class of "quadratic" mono-scalar-tensor 
theories used in {3, 4), where the (field-dependent) coupling strength a{if) between the scalar 
field and matter contains two parameters: a{ip) = uq + fiocp. Every pair (acA)) represents a 
specific scalar-tensor theory of gravity. As discovered in {3), for certain values of /3o, neutron 
stars can develop a significant scalarization, even for vanishingly small ao, if their mass exceeds 
a critical (/Jo-dependent) value. For this reason, this class of gravity theories is particularly well 
suited to demonstrate how the limit pO] ) probes a new gravity regime that has not been tested 
before (see Fig. 2b). The specific parameters and EOS in Fig. 2b have been chosen for demon- 
stration purposes. A change in the EOS, for instance, would lead to a modification in the details 
of the functional shape of a a, but would not change the overall picture. 

Constraints on the phase evolution of neutron-star mergers So far, the best constraints 
on dipolar gravitational wave damping in compact binaries come from the observations of the 
milhsecond pulsar PSR J1738-H0333, a lAV^f^M^ neutron star in a tight orbit {P^ ~ 8.5 h) 
with a spectroscopically resolved white-dwarf companion {14, 23). However, as discussed in 
detail above, such timing experiments are insensitive to strong-field effects that might only 
become relevant in the strong gravitational fields of high-mass neutron stars. Consequently, 
the dynamics of a merger of a 2 M© neutron star with a "canonical" neutron star or a black hole 
(BH) might have a significant contribution from dipolar GW damping, leading to a modification 
of the orbital dynamics that is incompatible with the sophisticated GR templates used to search 
for GWs with ground-based GW detectors, like LIGO and VIRGO, {34). With the results 



on PSR J0348-I-0432, in particular with the limit given in eq. (20), this question can finally 
be addressed in some details. For this purpose, we decompose equation ( [T5| ) into the 2.5pN 
contribution, that is matched by an appropriate GR template, and the 1 .5pN contribution, that 
drives the phase evolution away from the 2.5pN dynamics. Following {38, 39), we introduce the 
dimensionless orbital angular velocity 

u = Mrih = nMfcw , (21) 

where /gw denotes the frequency of the GW and 

^ OM (X.X,.)3/^ 



To leading order, one then finds 



96 

Mm = — (m"/^ + Sm^) , (23) 
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where 

5/ X,Xs f (a.-asf 



The observed GW cycles in a frequency band [/in, /out] can be computed as follows: 

fdt= if/f)df=- —du. 



Consequently, the difference between the 2.5pN dynamics and the 2.5pN + 1.5pN dynamics is 
given by 
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where we made no assumption about the size of the value for S. For the LIGO/VIRGO band 
Min «; Moutjj Fig. 5 gives AA^gw for the LIGO/VIRGO detectors, for which a typical bandwidth 
of 20 Hz to a few kHz was assumed, as a function of la^ - orol for two different systems, a 
2/1.25 Mq NS-NS system and a 2/10 Mq NS-BH system. Concerning the NS-BH systems, we 
considered the class of alternative gravity theories where BHs are practically identical to GR, 
and consequently used aB = 0. For instance, this is the case in scalar-tensor gravity theories with 
negligible time dependence of the asymptotic scalar field (39). For the NS-NS system an ex- 
treme case is represented by the assumption that only the massive neutron star has a significant 
scalar coupling strength a a, while the lighter companion behaves like a weakly self- gravitating 
body, meaning ob = ^o- Besides this, for the NS-NS system we have also performed calcula- 
tions using a hypothetical most conservative (maximal AA^gw) value for the effective coupling 
strength of the companion B, which is ub = 0. However, such an assumption seems unphysi- 
cal for a non-zero ao, where aB is expected to approach ofQ (and not 0) for less massive stars. 
With the limit obtained from PSR J0348-I-0432 we find a conservative upper limit for the dipolar 
phase offset of ~ 0.5 (NS-NS) and 0.04 (NS-BH) cycles, an amount that would not jeopardize 
the detection of the gravitational wave signal in the LIGO/VIRGO band (37). 



Formation via a common envelope and spiral-in phase Common-Envelope (CE) evolution 
(40,99) in X-ray binaries is initiated by dynamically unstable mass transfer, often as the result of 
a high mass-transfer rate and a large initial donor/accretor mass ratio, qi = M2/MpsR > 1 . If the 
CE is initiated while the donor star is still early in its main sequence stage (i.e. if Pb < 1 day), the 
outcome is expected to be a merger (100). It is generally believed that a binary can only survive 
the CE evolution, and thereby successfully eject the envelope of the donor star, if the binding 
energy of the envelope, ^bind^ is less than the released orbital energy from the in-spiral process, 

^For a detector that is sensitive up to a few kHz, the frequency /^^ is determined by the innermost circular 
orbit, which is ~ 1350 Hz for a 2/1.25Mq system and ~ 370 Hz for a 2/10 Mg system (see (35)). 
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AEorb(lOl). The orbital energy of PSRJ0348+0432 is: I^oaI = GM^s^M^^jla - S.SxlO^^erg. 
Hence, even if assuming in-spiral from infinity to the current orbital separation, the amount of 
liberated orbital energy from the CE phase cannot exceed this value. From calculations of 
EhmA of intermediate-mass stars (Table S3), we find that ^bind ^ ^E^rh during most of their 
evolutionary stages. Only if the donor star (i.e. the white dwarf progenitor) is an evolved giant 
is it possible to eject the envelope. However, in this case the core mass of such an evolved 
star, Mcorej is more massive than the observed white dwarf companion by at least a factor of 
2-3. (As argued in the main text, a reduction in white dwarf mass via evaporation from the 
pulsar wind seems to be ruled out for PSRJ0348-I-0432 and therefore cannot help circumvent 
this discrepancy.) 

From Fig. |7] we see that only a low-mass donor star with mass M2 < 2.2 M©, and not 
evolved beyond the terminal age main-sequence (TAMS), would leave behind A/core = ^^wd - 
0.17 Mq. In this case, it is clear that energy sources other than AEqA rnust contribute to expel 
the envelope (since in this case ^bind ^ A^orb)- Such an energy source could be the release 
of gravitational potential energy from material which accretes onto the neutron star during the 
CE. The amount of released energy per accreted unit mass is roughly tJJjm ~ GM/R ~ 2 x 
10^" erg g"^ Hence, assuming full absorption and 100% energy conversion of this released 
energy to eject the envelope, this would require accretion of ~ 4 x 10~^ Mq; a value which is 
not unrealistic given a timescale of the CE event of ~ 10^ yr with Eddington limited accretion 
(a few 10~*^Moyr-'). 

As a consequence of this relatively short CE phase, the currently observed mass of MpsR = 
2.01 M© should be close to the original mass of the neutron star after its formation in a type Ib/c 
supernova. According to recent studies by (102), neutron star birth masses of 2.0 Mq are indeed 
possible. As mentioned in the main text, however, having an initially massive neutron star 
would be a more serious problem for formation via a CE event with a < 2.2 Mq donor star. 
Such a high value of MpsR would lead to a value of qi close to unity, in which case the Roche- 
lobe overflow (RLO) is expected to be dynamically stable, thereby avoiding the formation of 
a CE (42, 43). The only solution to this problem would be that the neutron star was originally 
born with a more typical mass of ~ 1 .4 Mq, in which case qi would be sufficiently high to ensure 
formation of a CE. However, in that case one would have to accept the concept of hypercritical 
accretion (46, 103), allowing the neutron star to accrete a large amount of mass ~ 0.5 - 0.7 M© 
on a timescale of ~ 10^ yr. One could argue that PSRJ0348 -1-0432 would then be the best (and 
to our knowledge the only) candidate known in which hypercritical accretion might have been 
at work. 

To summarize, given the many issues discussed above we find that a CE formation channel 
is less favorable to explain PSRJ0348-I-0432 and we now proceed with investigating another 
solution, the LMXB formation channel. 

Formation via a converging low-mass X-ray binary As mentioned in the main text, a hand- 
ful of binary pulsars exist with values of Pb ^ 8 hr and Mwd « 0. 14 - 0. 18 M©, similar to those 
of PSR J0348-I-0432. These systems are tentatively thought to descend from low-mass X-ray 
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binaries (LMXBs) in which the binary suffered from loss of orbital angular momentum caused 
by magnetic braking (43, 47, 48). However, there remains a general problem for reproducing 
these pulsar binaries using current stellar evolution codes. A main issue is that converging 
LMXBs most often do not detach but keep evolving with continuous mass transfer to more and 
more compact systems with P^ <\\\x: and ultra-light donor masses M2 < 0.08 M©. In a few in- 
stances, where fine-tuning may lead to detachment and the right values of P^ and Mi, the donor 
star is typically too hydrogen rich to settle and cool as a compact He white dwarf [however, see 
sequence d in fig. 16 of {43) for an exception]. Our numerical studies are no exception from 
this general picture. 

Using the Langer stellar evolution code [e.g. {41, 104)] we have attempted to model the 
formation and evolution of the PSRJ0348 -1-0432 system. Here we present a solution where 
we have forced the donor star to detach its Roche lobe at P^ ~ 5hr, such that the system 
subsequently shrinks in size to its present value of PoA - 2.46 hr due to GW radiation within 
the estimated cooling age of the white dwarf (%d ^ 2 Gyr, depending on cooling models and 
assumed metallicity). To be more precise, the estimated %d is actually a lower limit on the 
timescale during which the detached system evolved via GW radiation since it takes 10^ - 10^ yr 
for the detached pre- white dwarf to settle on the final cooling track. This can be compensated 
for by choosing a slightly larger P^, at the ZAMS, which causes the system to detach from the 
LMXB in a somewhat wider orbit. 

In Fig. [8] (and see also Fig. 6 in the main text) we show an example of our LMXB calcula- 
tions. The model binary shown here consisted initially of a 1.75 M© neutron star and a 1.1 M© 
donor star with metallicity Z = 0.02, mixing length parameter, a = 2.0 and ZAMS orbital 
period, P\, = 2.55 days. The initial P^ depends on the modeling of magnetic braking. Here the 
value corresponds to onset of RLO at Pb - 0.65 days, shortly after the donor star ceased central 
hydrogen burning. Our high value of the initial neutron star mass is motivated from studies 
which show that the accretion efficiency in LMXBs must be rather small — even for systems 
which are expected to have accreted at sub-Eddington levels {14, 15,42). Hence, by adopting an 
accretion efficiency of only 30% we need an initial high-mass neutron star in order to reach the 
present mass of PSR J0348-I-0432. Note, that some neutron stars are indeed expected to have 
been born massive [for a discussion, see {41) and references therein]. The outcome of our cal- 
culations would possibly have been somewhat similar by assuming an accretion efficiency close 
to 100% and starting with MpsR = 1.3 M©. To model the loss of orbital angular momentum due 
to mass loss from the system, we adopted the isotropic re-emission model {105). 

Based on its proper motion and radial velocity measurements, PSR J0348-I-0432 has an esti- 
mated 3D space velocity of 56 ± 8 km s"^ with respect to the Solar System. From Monte Carlo 
simulations of its past motion through our Galaxy [following the method described in {14, 106)], 
we find that this velocity corresponds to a peculiar velocity with respect to the local standard of 
rest at every transition of the Galactic plane of 75 + 6 km s~^ . This result is rather independent 
of the applied Galactic model. From subsequent simulations of the dynamical effects of the 
supernova explosion, we find that a relatively small kick magnitude of w < 150 kms"' was 
imparted to the newborn neutron star, by probing a broad range of values of the pre-supernova 
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orbital period and the masses of the collapsing naked He-core and its companion star (the white 
dwarf progenitor). 

Spin evolution of PSR J0348+0432 A peculiarity of PSRJ0348+0432, compared to other 
recycled pulsars with similar P^, and Mwd> is its slow spin period, P = 39 ms and its high value 
of the spin period derivative, P = 2.41xl0~^^ss"\ cf. the unusual location of PSR J0348+0432 
in both the /'/'-diagram and the Corbet-diagram (Figs. S9, SIO). In particular, the Corbet dia- 
gram clearly displays the unique characteristics of PSR J0348-I-0432 with a small P^ and a large 
value of P. 

During the LMXB phase, a pulsar is generally expected to accrete much more mass and 
angular momentum than needed to be spun-up to a few milliseconds (49). In the same process, 
its B-field should have decayed significantly — typically to values < 10*^ G. However, for some 
reason the B-field (~ 2 x 10'' G) remained relatively high in PSR J0348-I-0432. In contrast, the 
other known binary radio pulsars with similar values of Pb < 8 hr and Mwd ~ 0.14 - 0.18 Mq 
(e.g. PSRs J0751-I-1807 and J1738-I-0333), besides from the many black-widow-like systems, 
have low B-fields and spin periods of a few milliseconds, as expected from current theories of 
LMXB evolution. 



In Fig. 11 we have plotted the past and the future evolution of PSR J0348-I-0432. In the 
upper panel is seen the evolution of /'b. In the lower panel is seen the spin evolution of the 
pulsar assuming different values of a constant braking index 2 < n < 5. If the estimated cooling 
age of ~ 2 Gyr is correct (and adding to this value a pre- white dwarf contraction phase between 
RLO detachment and settling on the final cooling track, yielding an assumed total age of about 
2 - 2.5 Gyr) we can estimate that PSR J0348-I-0432 was recycled with an initial spin period of 
about 10 - 20 ms. This relatively slow spin could be (partly) caused by enhanced braking of the 
spin rate, due to the high B-field of the pulsar, during the Roche-lobe decoupling phase when 
the progenitor of the white dwarf ceased its mass transfer (50). If the total post-LMXB age is 
~ 2.6 Gyr then the pulsar could, at first sight, have been recycled with an initial spin period 
of 1 ms for n > 3. However, calculations of the pulsar spin-up line (49) do not predict such a 
rapid spin for pulsars with high B-fields and which accreted with typical mass-accretion rates 
of M < 10"^ M^dd (evident from both theoretical modeling of the LMXB RLO and observations 
of LMXB luminosities (107)). 
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Right Ascension (J2000) 

Figure 1: Finding chart for the PSRJ0348+0432 system and the comparison star used in our 
analysis (see text), created from the archived SDSS g' image. 
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Figure 2: Left: White dwarf cooling age (measured from the onset of the core contraction) 
when the temperature reaches T^s = 10000 K as a function of the total hydrogen mass of the star. 
Each line depicts a different total mass (from 0.155 to 0.185 Mq). For each model, hydrogen 
burning through the pp-chain at the bottom of the stellar envelope cannot be initiated below 
a critical envelope mass limit. As a result the white dwarf cools in a few Myr. For models 
below ~ 0. 162 Mq a temperature of 10000 K cannot be reached regardless of the envelope size. 
Right: Finite-temperature mass-radius relation (for 10000 K) for models that have the minimum 
envelope mass required for hydrogen burning (red line). Over-plotted are the most-likely value 
and 1, 2 and 3cr constraints on the surface gravity for PSR J0348-I-0432 (solid, dashed and dotted 
black lines respectively). For masses below 0.162Mo the radius is an extrapolation from lower 
temperatures (in dashed red). 



37 



21-; .. ■• ■.. . 



r=20.61±0.08 



3.0^ 
2.5 
g 2.0 ■ 
^ 1.5- 
.> l.Oh^; 
iH 0.5 
ce 0.0 
-0.5 



0.90 0.92 0.94 0.96 

MJD - 2455958 



*«»'i>j>^!fej'-S>jf»)i5W>S^ 



0.0 



0.4 0.6 

Orbital Phase, ch 



Figure 3: Photometric (upper) and phase-folded (lower) light-curve of the white dwarf compan- 
ion to PSR J0348+0432 in u', g' and r'. 
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Figure 4: Pulse profiles for PSRJ0348-I-0432 obtained with the WAPP spectrometers at fre- 
quencies of 1170, 1410, 1510 and 1610 MHz. Two full cycles are displayed for clarity. Their 
(almost) perfect overlap indicates that there is little pulse profile evolution between 1 170 and 
1610 MHz. The 1410 MHz pulse profile is the template used to derive all TOAs. The TOAs cor- 
respond to integer phases in this plot, which mark the maximum of the fundamental harmonic. 
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Figure 5: Top: Polarization profile of PSR J0348+0432 obtained with the EfFelsberg Telescope. 
The upper panel shows total intensity (/, black), the linearly polarized intensity (L, red) and the 
circularly polarized intensity (V, blue). The lower panel shows the position angle of L measured 
at pulse longitudes where L exceeds 2cr measured from an ofF-pulse RMS. The red line shows 
the resulting fit of a Rotating Vector Model (RVM), which indicates an "outer line-of- sight" 
(see (82) for details). Bottom: Map of the RVM parameters a (the angle between the spin axis 
and magnetic axis) and ^ (the angle between the line of sight and the spin axis). The green 
region corresponds to combinations of a, [, for which the RVM provides a good description of 
the polarimetry of PSR J0348+0432. Based on the polarimetry alone we would have a large 
uncertainty regarding a and ^. However, if we assume that during the accretion episode that 
recycled the pulsar the spin axis of the pulsar ^s aligned with the orbital angular momentum 
(which has an angle / = 4(K2 + 0^6 to the line of sight) then a - 45°. The minimum angle 
between the magnetic axis and the line of sight is then given by yS = ^ - or = -4° 8. 
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Figure 6: Post-fit residuals from the GBT (red), Arecibo (gray) and EfFelsberg (blue) TOAs, 
obtained with the timing model presented in Table 1 . Top: Residuals versus time. No significant 
un-modeled trends can be found in the TOA residuals. Bottom: Post-fit residuals versus orbital 
phase, which for this very low-eccentricity system is measured from the ascending node (i.e., 
the mean anomaly is equal to the orbital longitude). No significant trends can be identified 
in the residuals; indicating that the orbital model can describe the orbital modulation of the 
TOAs correctly. No dispersive delays or unaccounted Shapiro delay signatures are detectable 
near orbital phase 0.25 (superior conjunction), nor artifacts caused by incorrect de-dispersion 
or folding of the data. 
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Figure 7: Stellar core mass at different evolutionary epochs as a function of zero-age main 
sequence (ZAMS) mass. Assuming A/core = ^^wd - 0.17 Mq (as observed in PSR J0348+0432) 
constrains the progenitor star ZAMS mass to be < 2.2 Mq and that its envelope was lost near the 
terminal-age main sequence (TAMS). All calculations were performed without convective core 
overshooting. Including this effect (for example, using Sow = 0.20) would lower the required 
donor mass even more. [Figure adapted from (41)]. 
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Figure 8: Formation of PSR J0348+0432 from a converging LMXB for the same model as 
shown in Fig. 6 in the main text. The plot shows how P\^ (red line) and the mass of the accreting 
neutron star (green line) evolved as a function of decreasing donor star mass (here assumed to 
be 1.1 Mq on the ZAMS). The RLO was initiated when P^ ^ 16 hr and detached when P\, ^ 
4.9 hr. In this model the initial mass of the neutron star was assumed to be 1 .75 M©, although it 
may have been significantly lower if the neutron star accreted with an efficiency close to 100%. 
The present location of PSR J0348+0432 is marked with a star. 
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Figure 9: A /'/'-diagram of the 111 known binary radio pulsars in the Galactic disk. The 
location of PSR J0348+0432 is marked with a black star in a region which is mainly dominated 
by slow spin and high B-field pulsars with massive white dwarf companions (marked with blue 
diamonds). The dashed lines of constant B-fields were calculated following (49) and assuming 
for simplicity MpsR = 1 .4 M© and sin a = = a»c = 1 . AH P values in this plot are intrinsic 
values obtained from kinematic corrections to the observed values. Data taken from the ATNF 
Pulsar Catalogue (108), in Oct. 2012. 
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Figure 10: A Piy- P (Corbet) diagram of the 63 known Galactic binary pulsars with a He white 
dwarf companion of mass Mwd > 0.14Mo. The unique location of PSR J0348+0432 is shown 
with a star. Another puzzling pulsar, PSR J 1744-3922 (45) marked with a green asterisk, is 
included in this plot. These two pulsars seem to share high B-field properties with the 6 pulsars 
in circles. Data taken from the ATNF Pulsar Catalogue (108), in Sep. 2012. 
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Figure 11: Orbital evolution (top) and spin evolution (bottom) of PSR J0348+0432 in the past 
and in the future. Different evolution tracks are plotted for different assumed values of a (con- 
stant) braking index, n. The grey shaded region marks the estimated white dwarf cooling age. 
The past is for a negative value of time, future is for a positive value of time. See text for a 
discussion. 
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Table 1 : Observations log and radial velocity measurements 

Notes: (1) refers to the barycentric mid-exposure time. (2) is the orbital phase, (p, using the ephemeris 
in Table 1. (3) is the comparison's velocity in respect to the Solar System Barycenter and (4) the raw 
barycentric velocities of the white dwarf companion to PSRJ0348+0432. 



-1^ 

0\ 



Target No. 


^K.. 


0^ 


slit 


exposure 

s 


rotation 
deg 


pos. angle 
deg 


(kms-i) 


4d 
(kms^i) 


PSRJ0348+0432 1 


55915.070159 


0.9742 




799.96 


-135.0 


-150.0 


-18.24 + 0.41 


-348.25 + 7.17 


2 
3 


55915.080011 
55915.099359 


0.0704 
0.2593 




799.96 
799.95 


-135.0 
-135.0 


-155.1 
-166.5 


-22.46 ± 0.39 
-29.37 + 0.40 


-345.54 ± 7.55 
-12.36 + 9.71 


4 


55915.109036 


0.3538 




799.96 


-135.0 


-172.8 


-25.21+0.39 


+ 159.76 + 9.54 


5 


55915.120877 


0.4694 




799.96 


-135.0 


+ 179.2 


-19.26 + 0.38 


+332.63 ± 8.37 


6 


55915.130692 


0.5652 




799.98 


-135.0 


+ 172.7 


-20.79 ± 0.38 


+295.18 + 8.08 


7 
8 


55915.140364 
55915.150194 


0.6596 
0.7556 




799.97 
799.98 


-135.0 
-135.0 


+ 166.4 
+ 160.4 


-19.54 + 0.38 
-19.41+0.38 


+ 149.72 + 8.77 
-45.85 ± 8.76 


9 


55915.217387 


0.4116 




799.96 


-134.8 


+ 132.1 


-8.28 ± 0.44 


+311.34+ 11.94 


10 


55915.227349 


0.5089 




849.96 


-134.8 


+ 129.7 


-9.74 ± 0.43 


+360.61 + 11.87 


11 


55915.237602 


0.6090 




849.97 


-134.8 


+ 127.3 


-5.13 + 0.44 


+278.11 + 11.03 


12 


55915.247859 


0.7092 




849.95 


-134.8 


+ 125.4 


-4.85 + 0.46 


+86.48 + 13.56 


13 
14 


55915.261379 
55915.271051 


0.8412 
0.9356 




799.96 
799.97 


-134.8 
-134.8 


+ 123.1 
+ 121.7 


-3.10 + 0.42 
-2.20 + 0.46 


-143.98 ± 13.42 
-297.87 ± 14.87 


15 


55915.280729 


0.0301 




799.96 


-134.8 


+ 120.4 


-0.34 ± 0.52 


-272.12+18.87 


16 


55915.290404 


0.1245 




799.97 


-134.8 


+ 119.5 


+ 1.56 + 0.61 


-188.66 + 33.02 


17 


55916.060700 


0.6452 




799.98 


-134.8 


-146.8 


-49.62 + 0.44 


+ 197.14+ 15.88 


18 
19 


55916.070535 
55916.080364 


0.7412 
0.8372 




799.97 
799.96 


-134.8 

-134.8 


-151.6 
-156.8 


-52.40 + 0.50 
-50.10 + 0.51 


-10.28 + 21.58 
-209.03 ± 20.37 


20 


55916.091598 


0.9469 




799.97 


-134.8 


-163.4 


-48.75 ± 0.49 


-339.79 ± 20.05 


21 


55916.101421 


0.0428 




799.96 


-134.8 


-169.6 


-43.97 ± 0.44 


-354.42 ± 15.82 


22 


55916.111221 


0.1384 




799.97 


-134.8 


-176.1 


-41.06 + 0.42 


-255.63 ± 13.45 


23 
24 


55916.122672 
55916.132782 


0.2502 
0.3490 




799.96 
849.96 


-134.8 
-134.8 


+ 176.6 
+ 169.6 


-44.69 ± 0.46 
-46.23 ± 0.46 


-15.94+ 18.30 
+ 153.76+19.31 


25 


55916.142895 


0.4477 




799.96 


-134.8 


+ 163.1 


-43.02 ± 0.46 


+332.62+18.44 


26 


55916.154536 


0.5613 




799.96 


-134.8 


+ 156.3 


-41.91 ±0.49 


+285.43 ± 20.94 


27 


55916.164373 


0.6574 




799.97 


-134.8 


+ 151.1 


-40.10 + 0.48 


+ 165.98+18.92 


28 


55916.174210 


0.7534 




799.97 


-134.8 


+ 146.4 


-39.74 ± 0.48 


-54.07 + 20.13 



Continued on Next Page. 



Table 1 - Continued 

Target No. MJDbar slit exposure rotation pos. angle vr vwd 

s deg deg (kms"') (kms"') 







29 


55916.229170 


0.2900 


1" 


799.96 


-134.8 


+ 128.5 


-12.85 ±0.44 


+89.37 ± 15.28 






30 


55916.238967 


0.3857 


1" 


799.97 


-134.8 


+ 126.5 


-9.24 ± 0.43 


+300.61 ± 17.14 






31 


55916.250488 


0.4982 


1" 


799.96 


-134.8 


+ 124.4 


-9.52 ± 0.43 


+350.95 ± 18.00 






32 


55916.260274 


0.5937 


1" 


799.97 


-134.8 


+ 122.8 


-15.17 ±0.55 


+331.23 ±22.73 






33 


55916.270071 


0.6894 


1" 


799.97 


-134.8 


+ 121.4 


-18.58 ±0.62 


+89.53 ± 33.20 






34 


55916.279869 


0.7850 


1" 


799.97 


-134.8 


+ 120.2 


-19.10 ±0.76 


-110.32 ±50.99 






35 


55915.181590 


0.0641 


2'.'5 


799.97 


-134.8 


+ 144.0 










36 


55916.190536 


0.9148 


2'.' 5 


799.97 


-134.8 


+ 139.4 








EG 21 


37 


55915.031970 




1" 


21.99 


0.0 


-25.8 










38 


55915.034856 




2'.'5 


21.99 


0.0 


-24.5 








HD 49798 


39 


55915.350264 




1" 


2.00 


0.0 


+72.6 










40 


55916.343776 




1" 


2.00 


0.0 


+71.3 










41 


55916.346515 




2'.'5 


2.01 


0.0 


+72.3 








LTT3218 


42 


55916.352238 




2'.'5 


22.01 


0.0 


+58.8 










43 


55916.357062 




2'.'5 


22.01 


0.0 


+62.4 






:^ 




44 


55916.369131 




1" 


35.00 


0.0 


+64.5 






"^ 


GD108 


45 


55916.360948 




2'.'5 


22.00 


0.0 


-176.4 










46 


55916.365961 




1" 


35.00 


0.0 


+ 179.9 







Table 2: Fractional binding energies of neutron stars (Data for Fig 4a). 



Neutron Star 


Mass (Mq) Reference Fractional Binding 


; Energy 


pulsars with white dwarf companions 


PSRJ0348+0432 
PSRJl 141-6545 
PSRJ1738+0333 


2.01 (this paper) 
1.27 {109) 
\A1 {14) 


-0.1446 
-0.0838 
-0.0993 


pulsars with neutron star companions 



PSR J0737-3039A 


1.338 


{22) 


PSR J0737-3039B 


1.249 


{22) 


PSRB1534+12 


1.333 


{110) 


. . . companion 


1.345 


{110) 


PSRB1913+16 


1.440 


{111) 


. . . companion 


1.389 


{111) 



-0.0890 
-0.0822 
-0.0887 
-0.0896 
-0.0969 
-0.0929 



Fractional binding energies of neutton stars in relativistic binaries, which are currently used in precision tests for 
gravity, and where the neutron star masses are determined with good (< few %) precision. The masses are taken 
from the given references. The specific numbers for the fractional binding energy are based on the 
equation-of-state ".20" of {27). A different equation-of-state gives different numbers, but does not change the fact 
that PSR J0348+0432 significantly exceeds the tested binding energy range. 
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